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Abstract 

A discontinuous area-preserving mapping derived from a sinusoidally- 
forced impacting system is studied. This system, the elastic impact 
oscihator, is very closely related to the accelerator models of particle 
physics such as the Fermi map. The discontinuity in the mapping 
is due to grazing which can have a surprisingly large effect upon the 
phase space. In particular, at the boundary of the stochastic sea, the 
discontinuity set and its images can act as a partial barrier which al- 
lows trajectories to move between chaotic and regular regions. The 
system at higher energies is also analysed and Moser's invariant curve 
theorem is used to find sufficient conditions for the existence of in- 
variant curves that bound the energy of the motion. Finally the be- 
haviour of the system under more general periodic forcing is briefly 
investigated. 
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1 Introduction 



We study a forced oscillator, such as a mass on a spring, moving in one 
dimension and repeatedly impacting against a fixed wall (Figure |l|). The 
motion between impacts is undamped and the impact is modelled as an 
instantaneous reversal of direction with a constant coefficient of restitution. 
For a linear spring with sinusoidal forcing this system, called the impact 
oscillator, is described after suitable rescaling by the following equations 

x + x = cos{ujt), x<a 
X — rx, X = a 

The inelastic impact oscillator (r < 1), was first studied |T|, 0, ^, ^ as a 
model for many important engineering phenomena including the rattling of 
engine components and the behaviour of structures in earthquakes. Further 
work has been done on such aspects as the bifurcation structure and the form 
of the chaotic attractors [|, H, |3, |, 3 . Here we consider the elastic case r = 1 
where no energy is lost at the impacts. For engineering applications r is often 
close to 1 and the elastic limit is an efficient way of exploring large areas of 
the phase space. It is also possible to exploit the time-reversal symmetry of 
the elastic case to locate the periodic orbits. However the elastic case has 
unusual behaviour of its own, and that is the subject of this paper. 

The motion of the impact oscillator consists of a series of impacts sepa- 
rated by smooth motion between the impacts. It is natural to describe these 
systems in terms of the impact map, Pi, relating the state of the system 
at one impact to the state at the next. Each impact is described by the 
(positive) velocity of the mass just before the impact and the phase (from 
to 27r) of the forcing cycle at the time of impact so that the phase space is 
a half- cylinder. Impacting systems that have been studied before include the 
accelerator models of particle physics |10|, |ll|, [12], |13|, [1^, and billiards, 



see 



17, IS, ISl and references therein. 



Although a mechanical impacting system is clearly discontinuous when 
regarded as a continuous-time dynamical system, the resulting impact map 
may be analytic and indeed most of the literature concerns impacting systems 
with impact maps that are either smooth or have discontinuities deliberately 
introduced into a higher derivative. The map Pj is discontinuous because 
of grazes, or zero-velocity impacts, see trajectory B of Figure ^ where the 
mass approaches the wall, touches it, and is then pulled away again. Nearby 
trajectories (trajectory C) will either (just) hit the wall and be deflected 
by it or miss the wall and go on to impact at a later time (trajectory A). 
Surprisingly, low velocity impacts strongly distort the phase space and so 
have a large effect upon the overall dynamics. This process is described in §13 



and is an important source of chaotic behaviour. Previous work [0, M, 20, 21 



has shown that grazing can introduce additional bifurcations such as the 
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sudden appearance or disappearance of periodic orbits and the interruption 
of period-doubling cascades. 



The phase space of Pi resembles that of the Fermi map [1^, 0, ^ in that 
the motion becomes more regular as the velocity of the impacts increases. 
Figure H shows a phase space plot of Pi for parameter values a = 0.1, a; = 2.5 
which has clearly separated regions of phase space giving rise to chaotic 
and regular behaviour. In particular, at low velocities the phase space is 
filled by a stochastic sea which, as the velocity increases, is interrupted by 
regular islands associated with elliptic fixed and periodic points. Above the 
stochastic sea the behaviour is more regular, with bounding invariant curves 
and elliptic islands. Note that the upper boundary of the stochastic sea is 
the lowest bounding invariant curve. 

However for slightly different values of a and u we observe a new phe- 
nomenon. Figure ^ shows a single trajectory for the parameter values (7 = 0, 
u = 2.8. The trajectory starts in the stochastic sea, exhibiting a typical 
sensitivity to initial data, but then, at irregular intervals, leaves the sea and 
moves along a smooth curve before reentering the stochastic sea. There is, 
for this example, less of a distinction between chaotic and regular motion as 
both can occur on the same trajectory. This very unusual behaviour is due 
to the discontinuity introduced by the grazing. The grazing discontinuity 
set and its iterates form structures at the upper boundary of the stochastic 
sea which act as a partial barrier allowing trajectories to move between the 
chaotic and the regular regions. Partial barriers also occur for smooth sys- 
tems although the mechanisms are quite different (see MacKay et al. [^). 
The existence of trajectories that can behave both regularly and chaotically 
is of great importance to the understanding and modelling of this system. In 
an engineering context, it is also important that such trajectories involve a 
large number of high- velocity impacts (with associated high wear-rates) and 
is all the more surprising since they arise directly from low velocity impacts. 

The outline of the paper is as follows. In §H we define the impact oscillator 
and the impact map and in §^ we analyse the effect of the grazing disconti- 
nuity on the nearby dynamics. In §^ we examine the upper boundary of the 
stochastic sea and show how interactions between the grazing discontinuity 
set and invariant curves can allow trajectories to alternate between stochastic 
and regular regions. §^ describes the relationship between the elastic impact 
oscillator and other impacting systems such as the accelerator models and 
billiards. Then in §^ motivated by the links with the accelerator models, we 
turn our attention to the behaviour of the system at high velocities where the 
impact map is smooth. We show that the existence of bounding invariant 
curves has an interesting dependence on the parameter values and briefly 
look at the case of more general periodic forcing. 
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2 The impact oscillator and the impact map 



The impact oscillator is a mass on a linear spring, moving in one dimension 
subject to a sinusoidal forcing, which impacts against a rigid obstacle referred 
to as the wall (Figure |^). The impact itself is modelled as an instantaneous 
process with coefficient of restitution r, < r < 1. By rescaling the system 
so that the amplitude of the forcing and the natural frequency of the spring 
are equal to 1 we obtain equation (|l|) where u is the forcing frequency and a, 
the clearance, is the position of the wall relative to the equilibrium position 
of the mass in the unforced system. Thus the system has three parameters 
r, a and u. In this paper r = 1. 

If the impacting condition is removed then the system is just the forced 
harmonic oscillator — a linear system whose behaviour is completely under- 
stood. But the impact oscillator is highly nonlinear and cannot be regarded 
as a small perturbation of a linear or integrable system. It is convenient to 
replace the time, t, by the phase, (p, which is defined by = {ut) mod 27r. 
So we have the state vector u = (0, where v is the velocity x, and 
the phase space Q = 2ttS^ x (— oo, a] x M. The vector u evolves under the 
discontinuous fiow $. 

A common approach to studying a system with periodic forcing is to 
define a Poincare surface = {(0, x, f ) : = 0o} (different values of 0o give 
different maps but the choice is unimportant). This defines the stroboscopic 
map, Ps, which maps the position and velocity at the phase 0o of the forcing 
cycle onto the position and velocity one cycle later. Note that S5 is always 
crossed transversely since time is increasing. But it is the impacts themselves 
that are of interest and the map Ps tells us very little about these — for 
example, during a single forcing cycle there may be no impacts or there may 
be many and determining which is the case would appear to be a very hard 
problem. 

Instead we study the mapping Pi which maps one impact onto the next 
one (|l|, D, H). Pi is an extremely useful tool for the analysis of the impact 
oscillator. It contains all of the interesting dynamical information, since the 
impacts are the source of the nonlinearity and the motion between them is 
just that of a simple harmonic oscillator. 

The map Pi is discontinuous due to grazing trajectories. These are trajec- 
tories which pass through a point (0, cr, 0), known as a graze. Nearby orbits 
either have a low velocity impact close to (0, a, 0) or miss the wall and impact 
at a later time and so Pi is discontinuous in the neighbourhood of grazes. 

We now examine the map Pi and its discontinuities in more detail. 
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2.1 The impact map 



We define the Poincare surface S/ = {(0, : x = a,v > 0}. Tlie result- 
ing map takes tlie pliase and velocity at one impact to the phase 
and velocity (0„+i,f„+i) of the next impact. Conventionally, the velocity is 
recorded just before the impact and we find it convenient to not regard grazes 
as impacts. Thus the recorded velocity is always positive. More formally, 
the impact map, Pj, is defined by 

Pi : S/ ^ S/ where S/ = 2nS'^ x and Pi : (0„, Vn) t— > Vn+i) 

The phase space of Pi is therefore an open half- cylinder. Technically, the 
impact map is not a Poincare map since the surface is not everywhere 
transverse to the flow because of grazing trajectories. 

The precise definition of Pi is actually more cumbersome than we have 
described here. It is straightforward to show that since the motion between 
impacts is recurrent one impact must lead to another |^ ^ and so the map 
is defined everywhere. However, if |(j| < 1, it is possible for the mass to stick 
to the wall until the acceleration becomes negative and it moves away. This 
behaviour occurs with measure zero and if it occurs more than once on a 
single orbit, the motion must be periodic (in contrast, for the inelastic case 
r < 1, the mass can stick to the wall via an infinite sequence of bounces 
which brings the mass to rest in finite time and for some parameter ranges 
this behaviour can be extremely important). 

The free motion (between impacts) is described by the first equation of 
([H) which is just the forced harmonic oscillator and can be solved exactly. If 
uj 1 and the time and velocity of an impact are given by to and vq then 
the position and velocity of the mass at time t are given by 

x{t;tQ,Vo) = (cr — 7 cos cuto) cos(t — to) 

+ (— "^0 + ^1 sin uto) sin(t — to) + 7 cos ut (2) 



v{t;tQ,Vo) = {—Vq + ivy sinuto) cos(t — to) 

— (a — 7 cos uto) sm(t — to) — ^^7 sin ut (3) 

where 7 = 1/(1 — cu^). This solution is only valid for to < t < ti where ti is 
the time of the next impact given by the first solution of the transcendental 
equation 

x{ti;to,vo) = a (4) 

To solve the system numerically we only have to find the time of the next 
impact using a root-finding scheme. This is both quick and accurate since 
no numerical integration is involved although there is always the possibility 
of not detecting extremely low-velocity impacts. 
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Implicit differentiation of (0), (|^) and (^) gives an expression for the 
Jacobian derivative of Pj (see 

DPi{(Po,vo) = 

{NoS-voC)/vi S/v, 

{NoNjv^ - vo)S - {No + N^vo/vi)C N^S/v^ - C 

where S = sin(ti — to), C = cos(ti — to) and Ni = cos{ujti) — a, i = 0,1. 
The Ni are the accelerations of the mass just before the impacts. Note 
that t appears instead of on the right-hand side of (j^). This is purely for 
convenience. 

It is immediate that 

\DP,\ = ^ 

as stated in In fact this is equivalent to showing that Pi preserves the 
measure v dip dv and we can simplify the determinant by introducing the 
coordinate change z = v^. Writing Pj as a map from (0, z) to (0, z) gives us 
the Jacobian derivative 

DPi(0o,^o) = 

/ (iVo5-v^C)/v/ir \ 

(6) 

V 2(7VoiVi - ^/I^i)S - 2(iVov/5r+ iVi^)C {N^S - ^iC)/^^ j 
which has determinant equal to 1 and so is an area-preserving mapping. 

The elements of DPi depend on both (0o,'yo) and (0i,fi) and so the 
Jacobian is defined implicitly. Also the elements of (j^) become unbounded 
as fi — s> 0, that is, as the next impact tends to a graze. Pj is therefore a 
map from the half-cylinder to itself which is everywhere smooth except on 
the 1-dimensional set S, the discontinuity set, where 

S = lira I (00, fo) ■ -P[(0o, ^o) = (0, e) for some G 27rS'^| 

The set S can be thought of as the pre-image of the line v = 0. The sets 
S are included in Figures ^ and ^ and a more complicated example is shown 
in Figure |[ S" is not, as one might expect, a smooth closed curve that spans 
the phase space but instead consists of one or more smooth curve segments. 
This is because for part of the line v = the acceleration is positive and 
the mass sticks to the wall. Also there are trajectories that have consecutive 
grazes and these cause the smooth segments of S to connect to one another 
in a complicated way as shown in Figure For a detailed description of 
S see [0. An important result that we shall need in §|| and §|] is that 5* is 
bounded above. This means that for sufficiently high velocities the impact 
map will be smooth. 
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The time-reversal symmetry of the system means that the phase space of 
Pi is invariant under the transformation 

-t, (j)^2n-(j) (7) 

and is symmetric about the midhne = vr. The forward image of the hne 
f = is therefore the image of S under the time reversal symmetry (|^). We 
call this line W. We define further iterates and preiterates of these sets by 

The map Pi behaves very differently on each side of S. On one side of the 
curve the trajectories have low velocity impacts and the phase space becomes 
very distorted as can be seen from the singular Jacobian for f = 0. On the 
other side, the trajectory just misses the wall and hits it at a later time and 
the Jacobian does not have large elements. 



3 The dynamics of the impact map near S 

To illustrate the nature of Pi close to the discontinuity set S we imagine a 
line segment I of initial conditions in phase space which transversely crosses 
S (Figure Let / and S meet at point B and call the resulting two 
sections of I, J"*" and I . The particle motions corresponding to the points 
A, B and C are those shown in Figure |^. As we move along /"•" from A to 
B the image curve Pi{I^) is traced out and moves towards the line W and 
meets it transversely. The endpoint Pi{B) lies on the line W but is not 
(generically) the image of the point B under the symmetry transformation. 
As we continue from B to C the next impact is now a low velocity impact. 
So the curve Pi{I~) grows out of the line v = and the line / would appear 
to be split in two. However, the second iterate of /~ rejoins Pi{I~^) meeting 
it at Pi{B) and Whiston ^ showed that P]{I~) meets W tangentially. The 
side of S that does not map directly to low velocity impacts is called the 
non-grazing side and the side that does is the grazing side. The line Pi'i^I') 
is locally stretched by a factor of of (9(e~2) where e is the distance from S. 
This stretching is described by a square-root singularity. Together the cutting 
and stretching have important implications for the dynamics. Because grazes 
tend to collapse trajectories onto PF, W and its iterates strongly influence 
the overall dynamics. The stretching also means that periodic orbits which 
include low velocity impacts, f <^ 1, are highly likely to be unstable and 
so lie in the stochastic sea as can be easily seen by examining the Floquet 
multipliers derived from (j^) or (j^). 

We end this section with a few final comments about the map Pi. The 
impact oscillator system is unusual in that the preferred Poincare section 
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permits tangential intersections with the flow. This is not because there is 
no available everywhere-transverse section (E5 for example) but because the 
impact map is a much more powerful tool for studying the impacts which 
are the source of chaos in the system. It is the tangential intersections with 
the Poincare surface that define the sets S and W, but it is the low velocity 
impacts that are responsible for the cutting and stretching that is observed 
near W Because a; = o" is both the Poincare surface and the impacting surface 
these effects are seen together. 



4 Chaotic motion and partial barriers 

We now study the main chaotic region, often known as the stochastic sea. 
The two phase space plots, Figures and H, show the stochastic sea lying 
between the line v = and the regular region which exists at high velocities. 
It is our intention to explain the novel behaviour seen in Figure ^ where 
trajectories move between the stochastic sea and the regular region lying 
above it. To do this we must look carefully at the interaction between the 
discontinuity set and regular curves at the upper boundary of the stochastic 
sea. 



4.1 The discontinuity set and the stochastic sea 

As we have already shown, the set S" is a very strong source of chaos in 
the system due to the arbitrarily large stretching close to a graze. From 
numerical experiments it appears to be impossible for an invariant curve to 



cross S without being destroyed. This agrees with previous studies ||TT], |T9 
which observed that invariant curves do not usually survive crossing lines of 
discontinuities in the first or second derivative and implies that S can only 
lie within chaotic regions. It is easy to obtain an upper bound on S from 
equations (|^) and (0) and also straightforward to show that 5* must touch 
the line v = (for all r, a and rational u there exists a smooth motion which 
never hits the wall but repeatedly grazes it). These facts strongly suggest 
(but do not prove) that the set S will lie in a single bounded chaotic region 
extending upwards from v = 0, namely the stochastic sea. This is supported 
by the numerical studies. 



4.2 The boundary of the stochastic sea 

First let us consider the case where S lies well away from the boundary of 
the sea. This is illustrated by Figure |] where S lies well inside the sea. 
The boundary is the lowest invariant curve that spans the cylindrical phase 



8 



space and since the map is smooth in the neighbourhood of this curve the 
boundary is the same as those found in smooth area-preserving maps, with 
regular curves above it and a sticky chaotic layer just below it. 

We now examine Figure ^. This shows a single trajectory, with initial 
condition lying within the stochastic sea and followed for 20000 impacts. In 
that time it makes 3 excursions (each containing many hundreds of impacts) 
onto regular curves surrounding two different period-5 orbits. There is no 
smooth boundary to the stochastic sea, but a partial barrier which allows 
trajectories to move between the stochastic sea and the elliptic islands where 
the trajectory moves slowly around regular curves until it reenters the sea. 
This is because for these parameter values S does not lie within the stochastic 
sea, but touches the boundary and intersects the islands associated with the 
two elliptic period-5 orbits. 

In order to understand this mechanism we first look at how a trajec- 
tory moves from an elliptic curve into the stochastic sea. Figure || shows a 
schematic blowup of the last few iterates of the motion along the elliptic curve 
that is intersected by S (so only every 5**^ iterate is shown). The trajectory 
moves slowly down the curve until it crosses S. The trajectory is now on the 
grazing side of S and so the next impact (which is not shown) will be a low 
velocity impact which did not occur on the previous cycle. The trajectory 
is now cut, stretched and reconnected as described in §^ so that it returns 
to the neighbourhood of the elliptic curve having been stretched tangentially 
along the line and then very rapidly disappears into the stochastic sea. 

The opposite mechanism, where the trajectory suddenly jumps out of the 
stochastic sea onto the regular curves is of course just the time reversal of 
the above mechanism. So the corresponding diagram is Figure ^ with the 
order of the points reversed, S replaced by W and replaced by 5*^. It 
is not necessary for the system to be time-reversible in order for trajectories 
to be able to leave the sea, the existence of such trajectories can be deduced 
from the area-preserving property of the map — if trajectories can enter a 
region then trajectories must also leave that region. 

We refer to the boundary in this partial barrier since the dis- 

continuity set and its iterates form structures which severely limit the rate of 
transport (of phase space area) across them but do not stop it entirely. Par- 
tial barriers for smooth systems (see [^) include cantori and the turnstiles 
associated with periodic orbits. The mechanism described here is completely 
different and, unlike these other cases, the motion on one side of the barrier 
is regular. Some sections of the partial barrier allow trajectories to move 
upwards while others allow trajectories to move down. A section of the bar- 
rier is enlarged in Figure For this piece of the barrier the upward and 
downward sections are separated by the invariant curve which just touches 
S and surrounds the higher of the period-5 elliptic orbits. 
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The lower of the two period-5 elhptic orbits (the one that appears to 
he in the stochastic sea) does not actually exist, even though part of the 
surrounding curves do. It corresponds to an unphysical motion — that is, 
one in which the mass moves through the wall instead of experiencing a low- 
velocity impact. The position of these unphysical orbits can easily be found 
by modifying the numerical code so that low- velocity impacts are ignored and 
the mass is allowed to move freely for a short time in the forbidden region 
X > a. 

The precise nature of the partial barrier depends upon the number, type 
and orientation of the regular curves and regions that are intersected by S 
and it is possible to find parameter values for which the pictures can become 
extremely complex, for example when S intersects higher order chains of 
Birkhoff elliptic points. 

There is however another simple example that is very important, shown 



in Figure |I0|, and that is when S intersects a region filled with bounding 
invariant curves rather than a elliptic region. At first glance Figure |1^ looks 
just like Figure ^ with a stochastic sea bordered by a smooth invariant curve. 
However 5* now touches the boundary and the difference lies in the position 
and the nature of this invariant curve. It is the lowest bounding curve not 
because it is on the point of breaking up, with regular curves above it and 
a sticky chaotic border below, but because it lies immediately above S and 
all the invariant curves below it have been destroyed. Let us now consider 
the region just below this bounding curve. Here the phase space consists of 
invariant curves interrupted by S. Once a trajectory lands on one of these 
curves (having been stretched along W hy a. low velocity impact) it will move 
regularly along that curve until it lands on that small part which lies on the 
grazing side of S. It then has a low velocity impact and moves back down 
into the sea. This process is essentially the same as that of Figure ^ — the 
only difference being that in Figure |TD| the regular curves interrupted by S 
are bounding curves rather than elliptic curves and it is only the latter which 
rise above 5* and produce startling pictures such as Figure ^ So, while the 
stochastic sea in Figure ^ looks like a perfectly ordinary chaotic region from 
a smooth system it does in fact display the same alternating regular/chaotic 
behaviour as Figure ^. 

The grazing mechanism provides an interesting new transition between 
the regular curves that exist at high velocities and the truly chaotic motion 
at very low velocities. If we replace the wall by a steep potential gradient 
then the resulting smooth system should also display similar behaviour. 

We end by noting that exactly the same interactions between the discon- 
tinuity set and regular curves that occur at the upper boundary can occur at 
the boundaries between the chaotic region and elliptic islands that lie within 
it. 
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5 Related impacting systems 



Impacting systems and discontinuous mappings have been studied before. 
To help put our work into perspective we briefly look at two closely related 
classes of impacting system which have been fundamental in the study of 
both smooth and discontinuous dynamical systems. 



5.1 Billiards 



Billiards are an important class of impacting system [O, 18, 19|. A billiard 



is a point mass moving in a bounded 2-dimensional region with boundary V. 
When the mass hits V it bounces away elastically, following the usual 'angle- 
of-incidence equals angle-of-refiection' law. A boundary component is called 
concave or dispersing if it curves away from the bounded region and convex 
or focussing if it bends inwards. A billiard is called convex if the bounded 
region is convex. Natural, area-preserving, coordinates to use for the impact 
map are the curvilinear distance, rj, along the boundary from some arbitrary 
point and s = sin a where a is the oriented angle between the normal to the 
boundary and the incoming trajectory (see Figure The mapping from one 
impact {riQ, sq) to the next (r/i, Si) has an imphcitly defined Jacobian [|18], [T9 



J 



d{r]i,si) ^ 
dir]o,so) 

( Cod — cos ao 

COStti 



Co cos «! + Ci COS ao — CoCid 



cos ao cos ai 



Cid — cos ai 
cosao 



where Co, Ci denote the curvature at the points of impact and d is the Eu- 
clidean distance between the points of impact. The form of this Jacobian in 
the limit cosai ^ 1, which corresponds to the mass hitting a non-convex 
piece of the boundary almost tangentially, is the same as that of the impact 
oscillator with area-preserving coordinates in the grazing limit 2:1 -C 1. In 
other words, a low velocity impact in the impact oscillator corresponds to a 
nearly tangential impact in a non-convex billiard which will have a disconti- 
nuity set at the pre-images of such impacts. In fact the cutting and stretching 
close to grazing trajectories described in ^ and Figure ^ is exactly the same 
for these non-convex billiards. 

The most frequently studied billiards are those with boundaries that are 
either every where-fo cussing or everywhere-dispersing. Billiards that have 
boundaries with focussing and dispersing components are likely to have both 
regular and chaotic regions and the mechanism described in §| will also be 
of importance to the study of these systems. 
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5.2 Accelerator models 



In 1949 Fermi proposed a mechanism for the acceleration of cosmic rays 
that involved collisions with magnetic field structures. Much work followed 
in which this process was modelled by particles impacting repeatedly (and 
elastically) against oscillating heavy objects. The idea also found a natu- 
ral application to particle accelerators where imperfections in the confining 
plasma ring could be modelled as periodic 'kicks' in a very similar manner. 
Because of these applications, and the simple form that the mappings take, 
they have become standard problems in Hamiltonian systems with two de- 



grees of freedom. Some work has also been done for inelastic systems p5| , |26 

Several different accelerator models have been studied. Two of the most 
important are the Fermi model and the Pustyl'nikov model. 



5.2.1 The Fermi models 



Ulam et al. |^ performed the first study of the Fermi mechanism. The 
model they used was a particle bouncing elastically with constant speed be- 
tween two walls — one fixed and one moving periodically. The impact maps 
arising from such systems are known as Fermi maps. They were primarily in- 
terested in the long-term behaviour of the particle and the maximum velocity 
that could be attained. Further studies can be found in [0, |1^, [1^, . 



Typically the grazing discontinuity set lies far below the lowest bounding 
invariant curve which is why the chaotic/regular trajectories described in 
§ ^ have not been observed for this model. A simplification that has been 
frequently made is to treat the position of the oscillating wall as being fixed 
while allowing its velocity to oscillate. In this way, the time of the next impact 
can be calculated explicitly for many wall motions and so the mapping itself 
can be written explicitly. This is a good approximation for high-velocity 
motions although at low velociies it precludes the possibility of grazing. For 
a sinusoidal wall motion the simplified mapping is 

Un+i = |M„ + sin(^„)| 

^pn+i = {ijJn + A/ Un+i) mod 27r 

where A is the system parameter. Note that the time between impacts is 
inversely proportional to the velocity. 

At low velocities there is a stochastic sea interspersed with elliptic islands. 
At higher velocities it has been proved (see |T^) that for wall velocities that 
are C^~^'^ there are invariant curves which cross the phase space and act as 
upper bounds to the velocity of the particle. If the wall motion is given by, for 
example, a saw-tooth function then these bounding curves disappear and the 
motion can become unbounded. Indeed much of the subsequent work on this 
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system was concerned with the existence of such curves and the smoothness 
conditions of KAM theory. 



5.2.2 The Pustyl'nikov map 

Pustyl'nikov [ITHI examined a different system consisting of a ball returning 



to an oscillating wall under the influence of gravity. Again the position of 
the wall can be regarded as fixed for high velocities and this gives rise to the 
following simplified mapping 

Un+l = \Un + Sm{lpn)\ /gx 
= {ipn + AUn+l) mod 271 

for sinusoidal forcing. Pustyl'nikov proved that even for analytic wall veloci- 
ties parameters can be chosen such that there is no bound to the velocity of 
the particle. The difference between this case and the Fermi maps is due to 
the time between impacts at high energies which now is proportional to the 
velocity. 



5.2.3 The impact oscillator 

In (|l|) the wall is stationary and the forcing is on the spring. However it is an 
easy exercise to show that, upon making the substitution y = x — •ycos^ut), 
equation (|I|) becomes 

y + y = 0, y < a - cos{ujt) 

y \—¥ —ry + {l + r)^ujsm{ujt), y = a — 'y cos{ujt) 

This is the equation of motion for an unforced mass moving on a linear 
spring and impacting against a wall whose position at time t is given by 
y = a — 7 cos (cut). Therefore, except for = 1, the forced oscillator problem 
defined by ([l|) is equivalent to a moving wall problem. This equivalence 
also holds for more general periodic forcing as long as there is no power at 
the frequency a; = 1. This relationship between the two systems is highly 
relevant since in physical situations the driving oscillations may act on the 
wall or on the mass/spring component. 

In the Fermi model the particle moves with constant speed between im- 
pacts with the vibrating wall. In the Pustyl'nikov model the particle moves 
freely under a constant acceleration (linear potential). The next 'natural' 
case to consider is a particle moving in a quadratic potential which, by the 
above result, is the impact oscillator (for u ^ 1). So the impact oscillator 
fits very naturally into the family of accelerator models. 

The impact oscillator is similar to the Fermi model with chaotic behaviour 
at low energies and regular behaviour at higher energies. But there are 
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important differences, especially at high velocities. These are also due to the 
time between impacts. For the Fermi model the time between impacts is 
0{l/v) for large velocities while for the Pustyl'nikov model it is 0{v). For 
the impact oscillator the time between impacts tends to tt as f increases as 
can be seen from equations (H) and (|^). So at high velocities the period 
of the oscillation is almost independent of the amplitude. This property is 
reminiscent of linear oscillators and for this reason we expect resonance effects 
to be important. This is indeed the case and the behaviour at high velocities 
is much more complicated than for the Fermi or Pustyl'nikov models, as we 
show in §^ 



5.3 Discontinuities in impacting systems 

We briefly discuss the different kinds of discontinuities that can appear, and 
have been studied, in impact maps. 

A billiard with a smooth boundary has a smooth impact map and in- 
troducing a discontinuity into the n^^ derivative of the curvature at some 
point of the boundary leads to an impact map with discontinuous (n — 1)**^ 
derivative. This situation has been much studied, especially with respect to 
proving ergodic and mixing properties and testing the smoothness conditions 
of KAM theory. Similarly, if for the simplified Fermi map, the wall velocity is 
smooth, then the impact map is also smooth and introducing discontinuities 
into the wall velocity or a higher derivative of the wall velocity results in a 
corresponding discontinuity in the impact map. 

When the map itself is discontinuous, nearby trajectories that straddle the 
discontinuity are separated. When the discontinuity is in the first derivative 
or higher, nearby trajectories remain close together. For examples of such 
discontinuous systems, see 0, [T^, n7\, ESI 



The grazing discontinuity is of a different form. Nearby trajectories that 
straddle S are not separated. Instead, one trajectory has an extra (low 
velocity) impact and the distance between the trajectories is stretched as 
was described in the previous section. It is perhaps helpful to think of a 
graze as being less severe than a discontinuity in the mapping but more 
severe than a discontinuity in the first derivative (if the extra low velocity 
impact is ignored then the impact map is Holder continuous with exponent 
|). Grazing discontinuities are a natural feature of impacting systems, much 
more natural than a discontinuity in the ra*^ derivative, but have been less 
thoroughly studied. 
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6 Behaviour of the impact oscillator at high 
velocities 



In this section we examine the impact oscillator at high velocities in the 
region above S where Pi is smooth. 

We showed in §H that the elastic impact oscillator fits very naturally 
into the family of accelerator models that includes the Fermi maps and the 
Pustyl'nikov maps. This link motivates us to ask the following question — 
under what conditions is it possible for the velocity of the mass to become 
unbounded? To answer this it is necessary to find conditions for the existence 
or non-existence of bounding invariant curves. 

Similar questions have also been asked of smooth systems. The bound- 
edness or otherwise of a particle moving in a one- dimensional smooth time- 



periodic potential has been investigated by several authors, see Norris |^ 
and references therein. Norris gave a sufficient condition for the existence of 
bounding curves for a system of the form 

x + g{x)=p{t) (11) 

where p{t) is a sufficiently smooth periodic function, g{x) oo as x —>■ ±00 
and either 

g{x)/x — > 00 as X ±00, or 

g{x)/x — » as X — * ±00 

The method of proof relies on the construction of a twist map and breaks 
down for functions g{x) which are 'too close' to the linear oscillator. The 
obvious example is the linear oscillator itself g{x) = x for which (|ll]) has un- 
bounded solutions whenever p{t) has a Fourier component with the resonant 
frequency 1. 

The impact oscillator retains some of the characteristics of the linear os- 
cillator which describes its motion between impacts, and for high velocity 
motions the time between impacts tends to a constant. Therefore, the fre- 
quency is approximately independent of the amplitude of the motion and we 
can expect there to be resonant forcing frequencies. 

We first consider the response of the impact oscillator to sinusoidal forcing 
and show that it depends not just on the forcing frequency but also on the 
position of the wall. We use the impact map Pi to examine (|l]) at high ener- 
gies and give sufficient conditions on the parameter values for the velocities 
to be bounded for all time. We also find parameter values for which there 
exist trajectories whose velocity can become unbounded. We then briefly 
look at the case of more general periodic forcing. 

We start by considering the approximate behaviour of the system at large 
velocities. 
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Lemma 6.1 For uj ^ 1 and v ^ max(l,u;7) the impact map has the follow- 
ing form 

Vl = Vo + f{(Po) + Oi-) 

cj), = (0o + a + ^ + O(^))mod27r (12) 

where a is a constant independent of{(j)Q,VQ). 

Proof We assume that Vq is large and, in particular, that the point (00? "^o) 
lies above the set S which is bounded. From equation (^, ti satisfies 

a = (cr — 7 cos uto) cos(ti — to) — ("^^o — ^7 sin uto) sin(ti — to) + 7 cos uti (13) 

and for large Vq the time between impacts is close to vr. So we let ti— to = 7r+(5 
where 6 is small and expanding (|13D in powers of l/i;o gives 

2cr - 7(cosa;to - cosu;(to + tt)) .^,1. 
d = + 0{—) 

where g{t) as defined in the statement of the Lemma is 

5'(t) = 2cr — 7(cosct;t + cosct;(t + vr)) (14) 

Substituting this value for ti into equation (^ we obtain the following ex- 
pression for Vl 

Vl = fo — ^7(sinu;to + sinu;(to + tt)) 

+ ^^[a - 7coscjto + (1 - 7) coscj(to + tt)] + 0{\) (15) 

We let 

/(t) = — ci;7(sinu;t + sinu;(t + tt)) (16) 
Replacing t by the phase in (|I4D and ([l^) we obtain 

Vl = Vo- cJ7(sin(0o) + sin(0o + a)) + C( — ) 

Vo 

01 = [00 + a + 

2a;a-^7(cos(0o) + cos(0o + a)) ^ ^1 ^vr (17) 

^0 v^ 

where a = ujtt mod 2tt. □ 
We now prove the following result. 

Theorem 6.2 If uj ^ 2n,n G and a ^ then the velocity of the impact 
oscillator is bounded for all time for all initial conditions. 
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Proof We prove this by showing that at high velocities the impact map 
is a perturbation of an integrable twist map of the form 

ri = ro 

^1 = ^o + /3(r) (18) 
where /?' 7^ 0. Then, if certain other conditions are satisfied, Moser's small 



twist theorem ^1[ guarantees the existence of infinitely many bounding in- 
variant curves. First we consider the case 7^ 1. 

We shall need the following Lemma. 

Lemma 6.3 Let na 7^ mod 2tt where n G N^. Then the functional equa- 
tions 

F{(f) + a) — F{(j)) = sin(n</) + na) + sin(n0) (19) 
G{(f) + a) — G{(l)) = cos{n(j) + na) + cos{n(f)) (20) 

are respectively solved by 

na 

F{(j)) = -cot(— )cos(n0) (21) 
G(0) = cot(-^)sin(n0) (22) 

Proof The result is easily established by using the Fourier transform. □ 

We now assume v ^ uj'-)coi{a/2) and make the following coordinate 
change 

w = u + u;7 cot(^) cos(0) 

ip = (f) ^cot(-)sin(0) (23) 

which, using the results of Lemma p.3| for n = 1, transforms (|l^ into 

wx = wo + 0{^) 

Wo 

V-i = ^o + a + ^ + 0(A) (24) 

Wo Wo 

Physically, this coordinate transformation corresponds to finding the approx- 
imate form of the bounding invariant curves, the existence of which we now 
prove. 



Now, following [31 1 we introduce a small parameter p via the transforma- 



tion w = p where < p < 1 and 1 < r < 2. Equation (^) becomes 

ri = ro + 0(p3) 

V^i = ^Po + a + ^ + Oip') (25) 
17 



This is a small perturbation of a mapping of the form ( ]T8| ) and we are 
now almost in a position to use Moser's small twist theorem. First we note 
that the map (^) is analytic and so easily satisfies the smoothness assump- 
tion of Moser's theorem. Secondly, if o" 7^ 0, then the twist condition /3' 7^ 
in equation (0) is also satisfied. It only remains to demonstrate the curve 
intersection property for (^), namely that any bounding curve C intersects 
its image C . This follows from the area-preserving property of Pj via the co- 
ordinate transformations we have made since. Even though at low velocities 
the phase space is cut and stretched by the set 5", the area under a bounding 
curve C which does not intersect 5" must be preserved. 

The result for uo ^ 1 now follows by a straightforward application of 
Moser's small twist theorem which guarantees the existence of invariant 
bounding curves for p^ ) for arbitrarily small p, which corresponds to ar- 
bitrarily high velocities of the original system. Therefore the velocity of the 
mass is bounded for all time. 

Finally we consider the case uj = 1. This is a special case because the 
equations of motion (j^) and (j^) do not apply. Instead the motion after an 
impact at (to, ^^o) is given by 

x{t) = crcos(t - to) - ('Wo + ^ sin(to)) sin(t - to) 

+ ^(t-to)sin(t) 
v{t) = -a sin(t - to) - ("^o + ^ sin(to)) cos(t - to) 

+ ^(t-to) cos(t) + ^ sin(t) 
Proceeding as before we find that for v ^ 1 the impact map has the form 

71 1 

Vl = vo- - cos(0o) + — ) 
2 vo 

1 TT 1 

= (/>o + 7r + — (2o- + -sin(0o)) + C(— )mod27r (26) 
Vo 2 Vq 

If we now consider Pj^ then (^) becomes 

V2 = Vo + 0{-) 

02 = 00 + — + 0(^) (27) 

and, for a 7^ 0, we can apply Moser's invariant curve theorem as before. This 
completes the proof of the thoerem. □ 
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6.1 Unbounded motion 



The resonant cases = 2n do in fact lead to unbounded motion for certain 
values of a. Since a = 0, (|17D becomes 

= t^o - 2to'7sin(0o) + C( — ) 

01 = 00 + — ((T-7cos(0o)) + mod27r (28) 

We prove the following result 

Theorem 6.4 For uo = 2n,n G N^, there are unbounded trajectories i/|cT| < 
|7|- 

Proof Making the substitution s = 1/v into (^) we get 

si = So + 2sou;7sin(0o) + C'(so) 

01 = 0o + 2socu(or-7cos(0o)) + O(s2) (29) 

We now define 0* by 0* = cos"-'^(cr/7) and 7sin(0*) < and let = 0j — 0*. 

Assuming that s,t 1 and calculating the first neglected term in (p9D 
we obtain the following 

ti = to(l + 2cu7Sosin(0*)) + O(s^to,soto) 

si = So(l + 2u;7Sosin(0*)) + O(s3,s2to) (30) 

By the definition of 0*, 2c(J7 sin(0*) < and so for sufficiently small s,t, 

l-Snl < l-Sn-il and < 

and it is clear that lim^^oo -s^ = which corresponds to a motion with 
unbounded velocity. □ 

We can summarise the above results as follows. For the generic case 
a ^ and u ^ 2n the motion is bounded at sufficiently high velocities. 
When a = the boundedness of the motion is marginal and depends on 
the higher order terms that were neglected in the above analysis. If the 
energy is in fact unbounded, then the rate of energy gain is extremely slow. 
The resonant frequencies are at uj = 2n where unbounded motion occurs for 
|cr| < I7I and numerical experiments indicate that this inequality is tight and 
for |a| > I7I the motion is bounded (in this case, the map (^9]) cannot be 
regarded as a small perturbation of an integrable map of the form (|TBp and 
Moser's theorem cannot be applied). 

Thus the existence of bounding invariant curves depends crucially on the 
position of the wall, a. For \cr\ sufficiently large (which means a ^ for 
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non-resonant u) the motion is bounded. In the mechanical system described 
by (|I]) there are two ways of changing the value of a. One way is of course 
simply to move the wall. The other is to add a constant term, a, to the 
forcing cos (cut). Then a simple coordinate change x ^ x — a recovers the 
initial impact oscillator system with the new value of a given hj a a ~ a. 
This leads us to the physically interesting, and rather counterintuitive, result 
that to suppress any unbounded behaviour it is enough just to add a constant 
force of sufficient magnitude — pointing either towards or away from the wall. 

6.2 General periodic forcing 

We can use the fact that the motion between impacts is linear to extend the 
method to the case of more general periodic forcing functions. We consider 
a forcing function p{t) with period 2tt/uj and Fourier components given by 

oo 

p(t) = ^ COS nut + g„ sin nut (31) 

n=l 

where, as described above, the constant term in the Fourier series has been 
taken as zero and absorbed into the position of the wall. We now prove the 
following result. 

Theorem 6.5 Let a ^ and p{t) G C^^^. Then for a set of irrational u 
having full measure the velocity of the impact oscillator is bounded for all 
time for all initial conditions. 

Proof In order to simplify the algebra we shall only consider an even 
periodic forcing function 

oo 

p{t) = ^Pn cos{nut) (32) 

ra=l 

since the terms of the sine expansion are dealt with in an identical fashion. 

The time of the next impact, ti, is given by the first solution of the 
equation 

oo 

cr = (cr - ^ 7„jo„ cos{nuto)) cos(ti - to)- 

n=l 

oo oo 

(^0 - ^^InPn sni{nuto)) sin(ti - to) + XI 7n.Pn cos{nuti) (33) 

n=l n=l 

where 7„ = 1/(1 - n'^u'^). \i p{t) e C^^' then the LHS of (|3|) is C^+' since 
n^n ~ '^/n as n — i> oo. So by the implicit function theorem tiitQ^vo) G C^^^. 
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Taking v ^ 1 and writing ti = to + it + 6 we obtain the following equation 
for S 



a 



(-cr + 7„p„ cos{nujto)) cos((5) + {vq - Y ^^InPn sm{nujto)) sm{6) 

n=l n=l 

oo 

+ InPn COs{nUj(to + IT + 6)) (34) 



n=l 



Proceeding as before, we get the following approximate mapping at suf- 
ficiently high velocities 



Vl 



oo -j^ 

f — ^i^7nPn(sinn<^o + sin(ri</)o + na)) + 0{ — 

^1 ^0 

00 + 

2auj - 'En=i t^7nPn(cosn0o + cos(n0o + na)) 



+ 



+ Oi-,) (35) 



where a = ujtt mod 27r and the mapping is C^^'^. Once again, for a 7^ 0, 
Lemma |6]^ suggests the correct coordinate change to reduce the mapping to 
a sufficiently small perturbation of an integrable twist map. This change is 



w 



Y nuJlnPn cot( — J COS[n(f 

n=l 



1 °° na 
— ^InPn cot(— ) sin(n</ 
^0 „=i ^ 



(36) 



and is only well-defined if both of the Fourier series in ( p6D converge. Ma- 
nipulating the coefficients of the first series slightly (if this series converges 
then so does the second one), we get 



E 

n=l 



ntU'JnPn COt( 



,na. 



< E 

n=l 
00 

= E 



nUJlnPn 



sm 



n=l 



2nU'JnPr, 



This last sum is the classic 'small-divisor' problem and for p G C^^^, con- 
vergence of the series is guaranteed if u is irrational and has a continued 
fraction expansion that satisfies certain conditions. The set of such u has 
full measure. Our smoothness condition on p{t) was chosen so that the per- 
turbation is C^'^'^ and we can now apply Moser's twist theorem exactly as for 
the sinusoidal case. □ 

So for a set of forcing frequencies that has full measure, the motion will 
be bounded for all C^"^"^ forcing functions (for a 7^ 0). It is not clear whether, 
for general periodic forcing, unbounded behaviour can still be removed by 
increasing the magnitude of the clearance. 
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From § |5.2.3| the impact oscillator with forcing 

oo 

p{t) = ^pncos(ncjt) 

n=l 

is equivalent to the moving wall problem with wall motion given by 

oo 

y(t) = 'jnPn cos{nujt) 



n=l 



So we have the immediate corollary that for the moving wall system the C^"*"*^ 
smoothness condition must be replaced by a C^"*"^ condition on y{t). This 
corresponds to a wall velocity which is C^'^'^. 
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Figure 1: The mechanical impact oscillator. 



Figure 2: Trajectories in the neighbourhood of a graze. 



Figure 3: Several trajectories of the impact oscillator showing the phase 
space of Pi for parameter values a — 0.1, lo — 2.5. The thick line is the 
discontinuity set. 



Figure 4: A single trajectory showing both regular and chaotic behaviour for 
(7 = 0, a; = 2.8. The thick line is the discontinuity set. 



Figure 5: The discontinuity set 5" for cr = 0, a; = 5.3. 



Figure 6: The dynamics close to S. 



Figure 7: Billiard coordinates 



Figure 8: Close-up of the partial barrier together with the last few iterates 
of a trajectory on an elliptic curve before it reenters the stochastic sea. 



Figure 9: An invariant curve, denoted by the solid line, separating pieces of 
the partial barrier that allow trajectories up and down. 
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Figure 10: The phase space for cr = 0, a; = 2.85. The discontinuity set is 
also shown. 
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